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We attempt the integration of the phenomenologically successful four zero texture 

of fermion mass matrices with the renormalizable SO(10) GUT. The resulting scenario 

en . 

is found to be highly predictive. Firstly, we examine the phenomenological implica- 
tions of a class of the lepton mass matrices with parallel texture structures and obtain 
interesting constraints on the parameters of the charged lepton and the neutrino mass 
matrices. We combine these phenomenological constraints with the constraints ob- 
tained from SO(10) GUT to reduce the number of the free parameters and to further 
constrain the allowed ranges of the free parameters. The solar /atmospheric mixing 
angles obtained in this analysis are in fairly good agreement with the data. 



1 Introduction 

The origin of the fermion masses and mixings along with the related problem of CP violation 
constitute a formidable challenge for elementary particle physics. Leaving apart extremely 
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small neutrino masses, even the charged fermion mass hierarchy ranges over at least five 
orders of magnitude. Since the fermion masses are derived from the Yukawa couplings, 
which are free parameters within the Standard Model (SM), these Yukawa couplings must 
span several orders of magnitude to accommodate the strongly hierarchical pattern of the 
fermion masses and mixings. However, the currently available data on the fermion masses 
and mixings is insufficient for an unambiguous reconstruction of the fermion mass matrices. 
To make matters worse, radiative corrections can obscure the underlying structures. Thus, 
the existing data cannot, without some additional assumptions, determine all the elements 
of the Yukawa coupling matrices for the quarks and the leptons. Some of these assumptions, 
invoked to restrict the form of the fermion mass matrices include the presence of texture ze- 
ros [1], requirement of zero determinant [2] and vanishing minors [3] to name just a few. The 
main motivation for invoking different mass matrix Ansatze is to relate the fermion masses 
and mixing angles in a testable manner which reduces the number of free parameters in the 
Yukawa sector. The recent evidence for non-zero neutrino masses and mixings leads to a 
further proliferation of free parameters in the Yukawa sector. In the absence of a significant 
breakthrough in the theoretical understanding of the fermion flavors, the phenomenological 
approaches are bound to play a crucial role in interpreting new experimental data on the 
quark and the lepton mixing. These approaches are expected to provide useful hints towards 
unravelling the dynamics of the fermion mass generation, CP violation and identification 
of possible underlying symmetries of the fermion flavors from which viable models of the 
fermion mass generation and flavor mixing could, hopefully, be constructed. 
The strong fermion mass hierarchy should be apparent in the fermion mass matrices them- 
selves with the contribution of smaller elements to physical masses and mixing angles ex- 
pected to be negligibly small. Thus, these elements can, effectively, be neglected and replaced 
by zeros: the so-called texture zeros. However, the current neutrino oscillation data is con- 
sistent only with a limited number of texture schemes p] . Specifically, the available neutrino 
oscillation data disallow all the neutrino mass matrices with three or more texture zeros [1] 
in the flavor basis. The texture zeros at different positions in the neutrino mass matrix, 
in particular, and the fermion mass matrices, in general, could be the consequence of some 
underlying symmetry. Such universal textures of the fermion mass matrices can be realized 
within the framework of Grand Unified Theories (GUTs). Though Grand Unification on its 
own does not shed any light on the flavor problem, the GUTs provide the optimal frame- 
work in which possible solutions to the flavor problem could be embedded. Grand unified 
models attempt to explain the masses and mixings in both the quark and the lepton sectors 
simultaneously. The textures for the mass matrices obtained in these models can either be 
assumed at the very outset or can be derived from the observed mixing matrix in the flavor 
basis. Alternatively, the textures of the mass matrix can be obtained by embedding some 
family symmetry within the chosen Grand Unification group. One particularly interesting 
class of models is that based upon the SO(10) Grand Unification group. There are two 



kinds of minimal models in this class: those based upon Higgs dimension 10, 126, 126 and 
possibly also 120 and/or 210 and those based upon 10, 16, 16 and 45 representations. The 
former choice, generally, has symmetric and/or antisymmetric texture mass matrices while 
the latter type generally imply lopsided mass matrices for the down type quarks and the 
charged leptons. In the present work, we consider the four zero texture [FTZ] Ansatz for 
the fermion mass matrices within the SO(10) GUT framework. Within this framework, we 
not only have a relation between the down type quark mass matrices (Md) and the charged 
lepton mass matrices (Mj) but also a relation between the up-type quark mass matrices (M u ) 
and the Dirac neutrino mass matrices (Mo)- 

Four zero texture Ansatz is especially important since it can successfully describe not only 
the quark but also the lepton sector including the charged lepton and the neutrino masses. 
Moreover, these are compatible with specific GUT models [4] and can be obtained from 
Abelian flavor symmetries j5]. Furthermore, these mass matrices can accommodate the 
present value of Sin2/3 [5J. In view of the phenomenological success of the FZT Ansatz, it 
would be interesting to examine it in the larger context of SO(10) GUT. The FZT Ansatz 
has been studied earlier [7] in the context of SO(10) GUT which, in general, leads to rela- 
tions between mass matrices. On the other hand texture zeros imply relations between the 
elements of the mass matrices. It is, therefore, interesting to examine FZT Ansatze within 
the framework of SO(10) GUT. In the earlier analysis [7j somehow, the constraints implied 
by texture zeros of the neutrino mass matrix have not been incorporated into the analysis 
and only SO(10) constraints on mass matrices have been used. In the present work, we 
have obtained the FZT constraints on the lepton mass matrices in Sec. 2. These constraints 
have been incorporated into the analysis of SO(10) GUT with four zero texture structure in 
Sec. 3. We have also presented the Takagi diagonalization of the general complex symmetric 
neutrino mass matrix in Sec. 3 without making any simplifying assumption for the phases. 
Conclusions have been summarized in Sec. 4. 



2 Four zero texture lepton mass matrices 

All information about the lepton masses and mixings is encoded in a hermitian charged 
lepton mass matrix M e and the neutrino mass matrix M v which is complex symmetric. 
Both these matrices have parallel texture structure i.e. zeros at (1,1) and (1,3) places [8]. 
In the present work, we consider a special case of FZT Ansatze in which the lepton mass 
matrices have parallel texture structure with zero entries at (1, 1) and (1, 3) places. The 
charged fermion mass matrices are assumed to be of the form 

/ a f \ 
M f = a} b f c f (1) 

V c) d f J 



where f — e,u and d for charged lepton, up-type quarks and down-type quarks, respectively. 
The elements a/ and c/ are complex with ay = |a/|e"^, c/ = \cf\e l ^ f . The elements bf 
and df are real because of hermiticity of M/. The unitary matrices V/ diagonalize the mass 
matrices M 
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diag(e i(t, f,l : e^f). 0/ is the 
diagonalization matrix for the real symmetric mass matrix Ml = PjMfPf. Using the 
invariants, TrMr, TrMr and DetM^ three elements (a/,&/,c/) of the total four elements 
can be written in terms of remaining one element df and the mass eigenvalues — ra/i , to/2 
and ra/3 of the fermion mass matrix M/ as 

b f = -m fl 
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The free parameter df should be in the range mf2 < df < m/3 for the elements |a/| and |c/| 
to be real. Using Eqn.(3) the elements of the diagonalization matrix, Of can be written in 
terms of the charged fermion masses and the charged fermion mass matrix element df as 
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The neutrino mass matrix 
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is complex symmetric for which all the four elements are in general complex and is diago 
nalized by a complex unitary matrix V^: 



M, 
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where M^ %ag = diag(— mi, 1712,1713) The diagonalization matrices V e , V u , V u and Vd are not 
physically observable but the products Uqkm = Kl^d [2]> an d Upmns — VeVv [10J are 
the physically observable quark and lepton mixing matrices, respectively. They can be 
parametrized in terms of the mixing angles and CP violating phases. Considering neutrinos 
to be the Majorana particles, we parametrize Upmns — U.P in terms of the three mixing 
angles (#12, #23 and #13) and the three CP violating phases (a, (3 and 5) as 
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iin 9ij and Cy = cos Oij for i,j—l, 2, 3. 
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The (1,1) and (1,3) entries of the neutrino mass matrix M v are assumed to be zero. Thus, 
we obtain two complex homogeneous linear equations for m\, 1712 and 7773 given by 



m\a + 7712b + m^c = 0, 

m\ad + m 2 6g + m 3 ch = 
where the complex coefficients a, b, c, d, g and h are given by 
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a — O e \\XJ e \ + O e i2U m i + OeisUtl, 

b = O ell U e2 + O el2 U m 2 + O el3 U t2 , 

C = O e nU e 3 + O e i2U m 3 + OeizUtz, 
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One can enumerate the number of parameters in Eqn.(ll). The nine parameters [three 
neutrino mixing angles (#12, #23, #13), three neutrino mass eigenvalues (mi, 777.2, ^3), two 
Majorana-type CP violating phases (a, (3) and Dirac-type CP violating phase, 5] come from 

and 



m T 



the neutrino sector and the four parameters [three charged lepton masses (m e , m 

free parameter d e ] come from the charged lepton sector, thus, totalling thirteen parameters. 

The three charged lepton masses are [TT] 



m e 



0.510998910MeV,m M = 105.658367MeV,m T = 1776.84MeV. 



(12) 
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Figure 1: Correlation plots for d e /m T , # 13 , M ee and the Majorana phases at 3a. 



The experimental constraints on the neutrino oscillation parameters with their best fit values 
and 3cr deviations are [121 
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In addition, there is an upper bound on the mixing angle #13 from the CHOOZ experiment (9° 2 < 
12.38). Recently T2K[13] and MINOS [14] experiments have given hints of a relatively large 

#13- 

For the simultaneous existence of two texture zeros at (1, 1) and (1, 3) positions in M u , the 
two values of mi given by 
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calculated from the mass ratios (^S ;^)> respectively, must be equal to within the errors of 
oscillation data. This constraint can be used to constrain the unknown parameters S, d e and 



the third mixing angle #13 which has not been measured experimentally as yet. We, thus, 
obtain not only the correlations [Fig. 1(a), Fig. 1(b)] amongst the unknown parameters [viz. 
the neutrino mixing angles, the CP violating phases and d e ] but also constrain the unknown 
parameter d e . We obtain two allowed regions for the free parameter d e . We call them 'low 
de and 'high d e ' solutions. The other two free parameters (0 e and ip e ) are not constrained by 
the presence of texture zeros because they factor out from Eqns.(8) and (9). We have, also, 
calculated the effective Majorana mass M ee (Fig. 1(c)) appearing in the neutrinoless double 
beta decay for the allowed parameter space 

M ee = \ mi U 2 el + m 2 U 2 e2 + m 3 U 2 3 \. (16) 

These constraints on the hitherto unknown parameter d e obtained from the FZT Ansatz will 
be used as an input in the next section where we attempt to constrain the SO(10) GUT. 



3 Four Zero Texture and SO(10) GUT 

In the next step we combine the FZT Ansatz with the constraints obtained from the SO(10) 
GUT. Such an analysis will require renormalization group (RG) running from weak scale to 
the GUT scale. However, it is known that the effects of RG running are negligible for mass 
matrices with normal hierarchy. Having obtained the texture zero constraints on the free 
parameter d e , we now further constrain it by imposing the following SO(10) GUT relations 
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where S, S' are the symmetric part coming from the 10 and 126 dimensional Higgs repre- 
sentations, respectively, and A is the antisymmetric part coming from the 120 dimensional 
Higgs representation. Here, 77, e, p, 7, 5', 5" and 5'" represent the relative coefficients of the 
vacuum expectation values (VEVs). From Eqn.(17), we find that eighteen parameters come 
from Sf + Af, f = u,d,e and four are the SO(10) coefficients 5', rj, e and 5", thus, totalling 
twenty two free parameters. Eqn.(17) implies the following relations 

4 V S u = (3 + V e)S d + (l- V e)S e , (18) 

5'A U = 5"A d = 5'5"A e . 



In the component form, the above set of SO(10) constraints (Eqn.(17)) can be written as 
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In the above, we have largely followed the notations and conventions of Ref . [7] which com- 
bines the SO(10) constraints in the form of Eqn.(19). We have 8 independent equations on 
22 free parameters appearing in the quark and the charged lepton sectors. However, using 
the experimental results on the quark/charged lepton masses and mixings which includes six 
quark masses, three CKM mixing angles, one CKM CP violating phase and three charged 
lepton masses, we are left with only one free parameter to be fixed. These eight independent 
equations can be used to relate d e , <p e and ip e with the corresponding parameters in the 
quark sector i.e. d u , d d , A0 = (p u — <p d and Aip = ip u — ip d (we note that only these phase 
differences are determined from the experimental data). The phases <p e and tp e are given by 

the equations 
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where the down type quark phases are 



tan i 



(1 - ^)c e 
a u sin A0 



tsaiip d 



ra d — a u cos A<p ' 

c u sin Aip 

rc d — c u cos Aip ' 



with r = 4- and 
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The parameter d e is constrained by the requirement that cos0 e and cos^e in Eqns. (20-21) 
are physical and 

r i «in in ., /-sin ih - 
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The parameters d u , d d , Acp and Aip can be constrained by the observed CKM matrix [6]. 
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Figure 2: Correlation plots for some of the SO(10) parameters for high d e 
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Figure 3: Correlation plots for phases from the charged lepton sector and d e /m T (higher 
region) . 

The best fit values realized for these parameters are given by 

(27) 

(28) 

(29) 

(30) 

In our numerical analysis, we have used the following best fit values for the quark and the 
charged lepton masses [7j: 

m u = lMMeV,m c = 302MeV,m t = 129GW, 

m d = 1.33MeV, m s = 26.5MeV, m b = l.OGeV, 

m e = 0.325MeV, m M = Q8.55MeV, m T = 1171.Q2MeV. 

The correlation plots between (f> e , ip e and d e /m T for the best fit values of d u , dd, A0 and Aif; 
have been plotted in Fig. (3). In the phenomenological analysis of the FZT Ansatz, presented 
in the previous section, we obtained two allowed regions for the free parameter d e compatible 
with the available data on lepton masses and mixings. Using these phenomenologically 
allowed ranges for d e along with the SO(10) GUT relations we constrain the d e — r plane. 
After the combined input of zero textures and SO (10) relations, the allowed ranges of d e are 
(0.307346m T < d e < 0.309016m T ) for the low d e and (0.935427m T < d e < 0.935991m T ) for 
high d e . So far we have succeeded in fitting the unknown parameters of the charged lepton 
sector for the best fit values of d u , dd, A0 and A^. 
The neutrino mass matrix from Type-I+II seesaw mechanism [15] is given by 

M V = M L - M D M^ l Ml (32) 
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which introduces three more unknown parameters p, 7 and 5'" (Eqn.(17)). The lepton mixing 
matrix is given by 

Upmns = V}V V (33) 

where V e and V v are the diagonalization matrices for the charged leptons and neutrinos. The 
diagonalization of the charged lepton mass matrix M e is similar to the quark sector (Eqn.(3)). 
For Majorana neutrinos, the complex symmetric neutrino mass matrix is diagonalized by a 
unitary matrix V v (Eqn.(6)). Takagi factorization [16] can diagonalize a general complex 
symmetric neutrino mass matrix without any simplifying assumption for the phases. For a 
complex symmetric neutrino mass matrix M v there exists a unitary matrix V v and a real 
non-negative diagonal matrix M„ 9 such that 

M v = V v Mf a9 Vj (34) 

where the columns of V v are an orthonormal set of eigenvectors for M V M},. The corresponding 
diagonal entries of M^ %a9 are non-negative square roots of the corresponding eigenvalues of 
M v Ml. The information about neutrino masses and mixings can be obtained from Eqn.(33) 
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In order to obtain consistent values of the lepton mixing angles we have to fine tune the 
unknown parameters p, 7 and 5'". For example, for d e = 0.93575m T , we obtain the following 
mixing angles 

#12 = 33.45°, (38) 

#23 = 47.92°, (39) 

#13 = 2.035°. (40) 

The unknown parameters p, 7 and 5'" are fine tuned to be 

p = 1.43 x 10~ 9 , 7 = 2 x 10 16 , 5'" = 10. (41) 

We obtain a small value for R v = A ^ 2 = 5.37 x 10 -4 , which is well below the experimentally 
allowed range. However, it must be noted that the mixing angles and mass-squared differ- 
ences have been calculated for the best fit values of the parameters Aifj, A</>, d u , da given in 
Eqns. (27-30). However, these parameters have wide ranges and one may obtain consistent 
values of #13 and R u when we consider the full ranges of these parameters. This has been 
illustrated explicitly in the Appendix by picking up some representative values of the param- 
eters from within their allowed ranges to calculate the neutrino mass squared differences and 
the mixing angles which are in agreement with the currently available data. Our objective in 
using the best fit values is only meant to illustrate the methodology presented in this work. 
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4 Conclusions 

We discussed the four zero texture quark/lepton mass matrices in the context of SO (10) 
GUT. The integration of these two scenarios reduces the number of free parameters to one 
when we incorporate the known data from the quark and the charged lepton sectors. In our 
analysis, firstly we obtained two regions of solutions for the unknown parameter d e from the 
phenomenological analysis for the four zero texture lepton mass matrices. Then we used the 
SO(10) relations to further constrain the allowed parameter space. We successfully fitted all 
the unknown parameters from four zero textures and SO(10) GUT within the framework of 
Type-I+II seesaw mechanism. The solar/atmospheric mixing angles obtained in this analysis 
are in fairly good agreement with the data. Of course, one must run down the neutrino mass 
matrix from the GUT scale to the electroweak scale and fit with the known oscillation data. 
However, it is well known that the effects of RG running are negligible for mass matrices 
with the normal hierarchy peculiar of the texture structure considered in this work. 
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5 Appendix 

In this appendix we give the numerical values of different parameters which satisfy the 
experimental data using the input parameters Aip, A</>, d u , d^ in there allowed experimental 
ranges[B]. Table-I lists the numerical values of different parameters. The neutrino mixing 
angles and mass squared differences for these values of parameters are 



6 12 = 32.26°, #23 = 


= 45.74°, 


013 = 12.34°. 


(42) 


Am 2 12 = 7.43 x lO^eV 2 , 


Am 2 32 


= 2.703 x 10~ 3 e1/ 2 . 
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